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ABSTRACT 


To  solve  an  Integer  linear  program,  we  identify 
particular  values  that  the  objective  function  can 
assume  on  fu^UJ.e  lattice  points,  inus,  we 
reduce  the  problem  of  finding  an  optimal  integer 
solution  of  n  dimensions  to  that  of  finding  a 
feasible  integer  solution  of  n-1  dimensions.  A 
Branch  and  Bound  Method  is  presented  to  solve  the 
latter  problem  for  the  0-1  case. 
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AM  ALGORITHM  FOR  INTEGER  LINEAR  PROGRAMMING 
BY  PARAMETRIC  MODIFICATION  OF  AN  ADDED  CONSTRAINT 

by 

Hajlme  Eto 

1.  INTRODUCTION 

We  define  the  standard  linear  programing  problem  to  be  that  of  finding 
a  vector  X  to 


(1.1) 


Maximise  z  ■  CX 


subject  to  the  constraints 

(1.2) 

(1.3) 


AX  <  B 

ai 

x  >  o 


where  A  la  mxn,X  is  n  x  1  ,  B  is  cxl  and  C  la  1  x  n  .  The 
integer  linear  programming  problem  is  a  linear  program  in  which  X  must 
also  satisfy 


(1.4) 


X  :  integral, 


i.e.,  each  component  of  X  must  be  an  integer.  We  will  call  a  vector  X 
satisfying  (1.1)  and  (1.2)  a  feasible  continuous  solution,  and  a  feasible  X 
tnat  also  maximizes  z  -  CX  an  optimal  continuous  solution.  Feasible  and 
optimal  Integer  solutions  are  similarly  defined  by  requiring  that  (1.3)  be 
satisfied. 

Rewriting  the  problem  (1.1),  (1.2)  and  (1.3)  in  the  following  manner: 


(2.1) 

(2.2) 

(2.3) 


Maximize  z  ■  CX 
Subject  to  AX  +  A*Y  -  B 

X,  Y  >  0 


3 

y 


and  applying  the  simplex  method  to  this  problem,  we  obtain  a  new  representation: 


0 


# 


2 

Maximize  z  ■  C°U  +  K 
Subject  to  A°U  +  IW  -  B° 

U,  W  >  0 

where  U  and  U  consist  of  the  same  variables  as  X  and  Y  but  in  a 
different  order.  We  call  the  components  of  Y  the  original  basic  variables 
and  the  components  of  W  the  current  basic  variables.  (W  and  Y  ,  or  X 
and  U  may  or  may  not  be  disjoint.)  In  an  optimal  representation,  C°  >  0 
and  B°  >  0  ,  and  an  optimal  continuous  solution  is  given  by  U  -  0  and 
W  -  B°  . 

General  Ideas  of  an  Algorithm  for  the  Integer  Programming  Problem 

Hereafter,  we  assume  an  optimal  continuous  solution  has  been  obtained, 
yielding  a  value  of  z  ■  zq  . 

Let  V  denote  the  set  of  all  the  variables  including  the  alack 
variables.  Let  us  denote  by  G  the  set  of  all  basic  variables.  We  may 
exclude  from  G  the  degenerate  variables,  i.e.,  the  variables  which  have 
zero  components  in  B  vector  in  the  current  tableau.  And  we  include  in  G 
the  dual  degenerate  variables,  i.e.,  the  variables  whose  shadow  prices  are 
zero  iu  the  current  tableau. 

Let  F  denote  the  set  of  all  the  variables  with  nonzero  coefficients 
in  the  objective  function.  In  many  integer  programming  problems,  most  of  the 
zero-one  integer  variables  have  zero  coefficients  in  the  objective  (see 
Discussion  2  below)  and  F  is  a  small  subset  of  V  .  If  a  variable  has  an 
upper  bound  less  than  one  in  a  range,  it  is  dropped  from  F  because  it 
must  take  on  value  zero  there. 

Therefore,  the  intersection  F  ft  G  is  a  relatively  small  subset  of  V  . 
Clearly,  F  O  G  varies  with  the  basis  exchange  because  G  does. 


We  can  assume  without  loss  of  generality  that  all  coefficients  of  the 
objective  are  integers  and  that  their  greatest  common  divisor  is  one. 
However,  the  greatest  common  divisor  d  of  all  coefficients  of  variables 
belonging  to  FOG  may  be  larger  than  one. 

One  of  our  basic  ideas  is  to  vt  the  feasible  convex  by  a  hyperplane 
parallel  to  the  objective. 

Experience  indicates  that  an  efficient  version  of  Gomory' s  method  of 
integer  forms  occurs  by  applying  Gomory  cuts  to  the  value  of  the  objective 
function  first.  (3)  (1)  (2).  In  the  Gomory  cutting  method  the  value  of  the 
functional  is  limited  to  be  integral.  But  from  a  well-known  theorem 
(Lemma  1  below)  in  number  theory,  the  value  of  the  functional  is  confined  to 
be  multiples  of  d  as  defined  above. 

After  a  candidate  for  a  permissible  value  of  the  functional  is 
identified,  we  search  for  a  lattice  point  at  which  that  value  is  attuined. 

In  most  integer  programs,  the  upper  bounds  of  unknowns  are  so  low  that  each 
of  them  can  take  only  a  few  values.  Therefore,  if  we  succeed  in  contracting 
the  permissible  range  of  values  of  the  variables  belonging  to  FOG 
sufficiently,  it  may  be  practicable  simply  to  enumerate  the  permissible 
combinations  of  integer  values  of  these  variables. 

The  other  basic  idea  is  to  irtroduce  the  dual  degenerate  variables  into 
the  basis  one  at  a  time  in  order  to  determine  a  more  restricted  range  of 
values  of  the  variables. 


■•*.  • 


“Vft*  *r**'  ii***. 
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2.  THEORETICAL  FOUNDATION  OF  THE  ALGORITHM 
Definition  1. 

Let  F  be  the  set  of  the  variables  which  may  change  the  value  of 
the  objective  function;  i.e.,  F  -  {x^  |  c^  +  0  and  u^  /  1} 
where  u^  ^  Xj  . 

Definition  2. 


Let  Gp  be  the  set  of  all  the  variables  whose  shadow  prices  are 
zero  at  an  extreme  point  P  whether  they  may  or  may  not  be 
basic  in  an  ordinary  sense  of  the  word  "basic."  If  a  degeneracy 
occurs  at  point  P  ,  all  of  the  degenerate  variables  can  he 
excluded  from  Gp  ;  i.e., 

Gp  ■  (Xj  |  Cj  ■  0  and  (the  associated  +  0  if  is  basic)} 


Definition  3. 


Let  dp  be  the  greatest  common  divisor  (G.C.D.)  of  these 
for  which  Xj  e  F  0  Gp  ;  i.e., 

dp  -  G.C.D.  (Cj  |  (Cj  i  0)  and  (Cj  -  0  at  a  point  P)  and  (u^  t  1)) 


Lemma  1 . 


If  c.  and  x^  are  integers,  then  \  c.x. 

3  j  3  J 


k  d  where  k  is  an 


integer  and  d  is  the  G.C.D.  of  Cj  .  The  proof  can  be  found  in  (4). 

In  what  follows,  some  preparation  is  made  for  Theorem  1  associated  with 
Figure  1. 

Let  Qq  be  a  point  inside  a  convex  polyhedron.  Let  H  be  a  hyper¬ 
plane  through  Qq  with  the  same  slope  as  the  objective  function  and  let  H 

•••»  Qn  • 


intersect  with  the  polyhedron  at  points  Q1,Q2» 


(Possibly 


Tifci'TT’i  'Pfv^'kflwKyHBH^. 


Qt  »  Qj  for  i  j*  j  if  more  than  two  constraints  meet  on  the  same  point.) 

Let  us  consider  the  smallest  subpolyhedron  containing  Qq*Q^»  •••* 

Let  ...,  Vy  denote  extreme  points  of  the  subpolyhedron. 

We  can  consider  another  parallel  hyperplane  H  that  intersects  the 


polyhedron  on  the  same  edges  as  H  does.  Let  Q.  ,Q-,  ...,  Q  denote  the 

1  /  N 


intersecting  points  so  that  and  are  on  the  same  edge  for  each  1 


Let  us  take  a  point  Qq  on  H  . 


Theorem  1.  (Theorem  of  a  Possible  Value  of  the  Objective  Function  to  a 
Feasible  Lattice  Point) 


If  there  is  any  feasible  lattice  point  L  in  the  subpolyhedron 
VjV2  ...  Vy  ,  then  the  value  of  the  objective  function  at  L  is  k  d^  , 
where  k  is  an  integer  and  d^  is  the  r.C.D.  at  the  ,-oint  (see 

Definition  3). 


Proof ; 

Let  Qj  and  be  on  the  edge  VjV^(j  +  k)  and  be  the 

basic  variables  at  Vj  ,  respectively.  and  can  be  expressed 

uniquely  as  a  linear  combination  of  variables  in  the  union  of  sets 
(Bj  U  B^)  .  A  similar  linear  combination  which  is  unique  holds  for  all  1  . 
Since  Qq  can  be  expressed  uniquely  as  a  linear  combination  of 
^1*^2*  Qfl-  »  Qq  can  uniquely  expressed  as  a  linear  combination  of 

(BXU  B2  U  ...  U  By)  . 

If  we  cut  the  polyhedron  by  a  hyperplane  H  parallel  to  the  objective 
function,  then  all  points  on  H  are  optimal  feasible  (not  necessarily 
Integer)  points  having  the  same  value  of  the  objective  function. 

Let  us  consider  an  extreme  point  of  the  cutting  polyhedron  on  H  ,  say 
Qx  which  is  assumed  to  be  on  the  edge  V^V2  .  By  repeating  the  process  of 
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appropriate  basis  exchange,  we  get  all  such  Q^Q^,  •••»  which  have  the 
same  value  of  the  objective  function.  This  fact  implies  that  the  variable 
introduced  into  the  basis  at  each  time  has  zero  shadow  price;  in  other  words, 
starting  at  all  the  variables  introduced  into  the  basis  belong  to  the 

se-  -  (ti^  U  I^)  •  Furthermore,  no  other  variables  belong  to  this  set, 
because  otherwise  H  would  have  other  extreme  points  having  different  values 
of  the  objective.  Hence  *  (B^  U  B2)  +  S^  where  is  the  set  of  all 

the  variables  whose  shadow  prices  are  zero  at  .  The  same  argument  holds 

for  H  and  and  leads  us  to  the  conclusion  that  G^  *  (Bj  U  B^)  +  ■ 

(B1  U  B2)  +  Sl  «  G:  . 

Similar  results  can  be  'btained  for  and  for  all  i  ;  e.g., 

suppose  G2  -  (B^  U  B^)  +  S2  *  ^2  *  ^3  U  ^4^  +  ^2  where  ^2  is  °n  e<*8e 

.  (B^  U  B^)  ,  the  basis  on  Q2  ,  is  obtained  from  (B^  U  B2>  uniquely 

by  introducing  a  variable  from  the  set  into  the  basis  and  deleting  a 

variable  from  the  set  (B^  U  B2>  whose  shadow  price  turns  zero;  in  other 

words,  a  variable  moves  from  to  (B^  U  B^)  and  another  variable  moves 

from  (B^  U  B^)  to  S2  while  the  others  remain  unchanged.  Thus  G]^  *  G2  • 
Furthermore,  we  see  tnat  Gi  "  "  G2  “  "2  8*nce  G^  «  G^  ,  anc*  ^2  "  ^2  * 

The  same  argument  leads  us  to  the  conclusion  that  G^  ■  G2  ■  . . .  ■ 

■  G^  «  G2  *  ...  ■  G^  .  Since  and  Qq  are  expressed  uniquely  as 

linear  combinations  of  Q^,  ...,  and  . QN  ,  respectively,  Gq 

and  Gq  are  the  same. 

Clearly,  since  tht.  variable  of  the  set  F  can  change  the  value  of  the 
objective  function,  we  can  drop  the  variable  which  do  not  belong  to  set  F 
from  consideration.  Thus,  the  theorem  is  proved  with  help  of  Lemma  1. 

Q.E.D. 


y  w* 
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Corollary  to  Theorem  1.  (Criterion  of  an  Nonexistence  of  a  Lattice  Point) 

If  t  0  (mod.  d^)  for  some  1  ,  then  no  lattice  point  can  exist 
in  the  subpolyhedron  specified  by  G  where  d^  -  G.C.D. 

(a^  |  Xj  e  O  G)  and  ■  {x^  |  a^  i  0  and  {  1}  in  (2.2). 


Proof; 


r  (i) 

From  Theorem  1,  I  . x  »  k  d  for  k  integer,  while 

4  I  „  larv  J 


j  |Xj  eFHG 


a . . x ,  ■  y  a.  ,x.  ■  b.  in  the  subpolyhedron  specified  by  G  .  Hence 
lj  J  jlxjCFnc  ij  j  1 

b^  ■  k  d^  for  k  Integer,  thus  b^  =  0  (mod.  d^)  for  all  i  if  any 


J 


lattice  point  exists  in  the  subpolyhedron. 


Q.E.D. 


Theorem  2.  (Theorem  of  Enumerability) 

If  there  is  any  lattice  point  in  the  subpolyhedron,  then  it  is 
necessarily  on  a  hyperplane  H  (or  one  of  the  hyperplanes  if  there  is  more 
than  one  such  H  in  the  same  subpolyhedron)  whose  value  of  the  objective 
function  is  k  d  (for  an  integer  k  ). 

Proof:  (See  Figure  2.) 

Let  us  assume  on  the  contrary  that  a  lattice  point  L  is  on  a 
different  and  parallel  hyperplane  whose  objective  function  value  is  not 
k  d  for  any  integer  k  .  Here  we  assume  without  loss  of  generality  that  an 
intersecting  point  H  of  H  with  an  edge  is  on  the  same  edge  as  .  M 
is  not  necessarily  a  lattice  point. 

Since  H  has  the  same  slope  as  the  objective  function,  every  point  on 
j]  has  the  same  value  of  the  objective  function.  Thus  the  value  of  the 
objective  function  on  M  rnd  the  value  of  the  objective  function  on  L  are 
equal  to  an  integer  value  v  .  The  basis  which  expresses  also 
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expresses  M  .  Thus  we  can  see 

(Gx  at  Q:)  -  (Gm  at  M)  . 

This  implies  v  is  equal  to  k  d  .  Hence  M  is  on  one  of  the  hyperplanes 
H  having  k  d  as  the  value  of  the  objective  function,  and  so  is  L  .  This 
co.itradicts  the  assumption. 

Q.E.D. 

Theorem  3.  (Theorem  of  the  Criteria  for  Entering  Another  Subpolyhedron  and 
Ending  the  Algorithm  Unsuccessfully) 

If  is  an  infeasible  point  before  a  basis  exchange  and  the 

algorithm  moves  to  a  feasible  point  Q^’  after  an  appropriate  basis 
exchange,  then  H  given  by  the  added  constraint  is  in  another  subpolyhedron. 
If  feasibility  can  no  longer  be  recovered  by  any  basis  exchange,  then  the 
added  constraint  is  out  of  the  feasible  domain. 

Proof  follows  from  Figure  3.  H  is  out  of  a  subpolyhedron  and  H  is 
out  of  the  feasible  domain.  Theorem  3  justifies  Step  3  and  Step  9. 
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3.  DESCRIPTION  OF  THE  ALGORITHM 
Step  1. 

Solve  the  standard  Linear  Programming  Problem  without  Integer 
requirement  (1.4)  and  obtain  zq  ,  the  optimal  value  of  the 
objective. 

Let  Pq  stand  for  the  optimal  extreme  point  of  the  continuous 
problem. 

Examine  whether  the  solution  is  integral  or  not. 

If  integral,  the  problem  is  solved. 

If  otherwise,  go  to  Step  2. 

Step  2. 

Define  d^  ■  G.C.D.  (a^  |  e  O  G)  for  all  i  for  a 
given  G  where  *  (x^  |  a^  +  0  and  u^  /  1}  . 

Examine  whether  b^  =  0  (mod.  d^)  for  all  i  or  not. 

If  not,  then  go  to  Step  9. 

Otherwise,  go  to  Step  3. 

Step  3. 

Add  to  (1.1),  (1.2),  (1.3)  and  (1.4)  the  following  constraint 

r  : 
o 

C  X  <  k  d 
-  o  o 

where  dQ  is  the  greatest  common  divisor  of  all  the  variables 


belonging  to  F  f)  G  . 
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Step  4. 

Examine  whether  or  not  the  solution  subject  to  the  augmented 
constraints  In  Step  3  Is  feasible  without  exchanging  the  basis. 
If  feasible,  go  to  Step  5. 

If  otherwise,  go  to  Step  9. 


Step  5. 


Obtain  the  Intersecting  points  of  the 

added  constraint  rQ  with  the  edge  of  the  convex  hull  by 
repeating  the  basis  exchange.  (One  of  the  points,  say 

=  (x^^.x^ . xi^  *  •••»  *N^)  *  is  already  obtained 

in  Step  3.)  In  more  detail,  exchange  an  appropriate  basic 
variable  with  one  of  the  dual  degenerate  variables  (i.e.,  a 
variable  whose  shadow  price  is  zero  in  the  current  tableau) 
and  get  another  intersecting  point  with  the  same  value  of  the 
objective. 

The  number  of  dual  degenerate  variables  is  equal  to  n-1  (or 
less  than  n-1  if  more  than  two  constraints  meet  at  the  same 
point),  where  n  is  the  number  of  the  components  of  X 
excluding  the  slack  variables. 

If  an  infeasibility  occurs  during  a  basis  exchange  (i.e.,  if 
one  of  . is  out  of  the  feasible  domain). 


go  to  Step  9. 


If  max  xj^  <  1  ,  then  define  F  ,  F  ,  d  and  dv~'  .  If 
j  (i) 

b^  ^  0  (mod.  new  d  )  for  any  i  or  the  current  value  of 
objective  is  not  k'd  for  new  d  and  any  integer  k'  ,  then 
go  to  Step  8  after  resetting  old  F  ,  F^  ,  d  and  d^  again. 
Otherwise,  go  to  Step  6. 


,(D 


0 
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Step  6. 


Obtain  a  lattice  point  (or  lattice  points)  that  can  be  expressed 

as  a  -?x  combination  of  the  points  ...»  Q.^  . 

Th,  s  that  the  x^  belonging  to  FOG  must  satisfy 

mx*  (x^.x*2* . xf>)  >  xt  >  min  (x^.xj2*.  -<2> 

,<«)  for  ail  Is  Flit  where  Q*0>  -  (x*J),x*J) 

x<J> 

m+n-1/ 


Step  7. 

Examine  whether  a  lattice  point  is  obtained  or  not. 
If  obtained,  the  algorithm  comes  to  an  end. 

If  otherwise,  go  to  Step  8. 


Step  8. 

Change  the  added  constraint  parametrically  to  r  ' 

C  X  <  (k  -  l)d 
■  o  o 

where  C  ,  X  are  the  same  as  in  Step  3  above.  Solve  the  problem 
subject  to  the  modified  constraints.  Go  to  Step  4. 


Step  9.  (From  Step  4) 

Identify  the  next  FOG. 
Go  to  Step  2. 


The  procedure  below  is  similar  to  the  above  case  if  we  substitute 

subscript  1  for  subscript  0  and  if  we  replace  "add  a  constraint  rQ 

C  X  <  k  d  "  with  "modify  the  added  constraint  r  into  r,  such  that 
=  oo  o  1 

r^  :  C  X  <  k^d^  •"  no  extreme  point  can  be  found  anymore,  then  the 
algorithm  terminates  unsuccessfully. 


,*..**«*sf3r»* 
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Remark  to  Step  2  and  Stei 


A  question  arises  how  to  know  In  advance  which  variables  belong  to 

F  H  G  after  the  addition  of  a  constraint  r  .  They  do  not  remain  the  same 

o 

as  on  Pg  . 

In  order  to  identify  the  variables  belonging  to  the  current  set  of 
FOG,  add  for  trial  a  constraint 


C  X  <  z  ' 

■  o 

where  z  '  ■  [z  ]  if  z  is  not  an  integer,  or  z  '  ■  z  -  1  if  z  is  an 
O  O  O  oo  o 

integer.  In  other  words,  z^'  is  the  maximal  integer  less  than  zq  . 

We  are  in  position  to  identify  the  variables  in  question  by  solving  the 
problem  under  the  added  trial  constraint.  Ti.e  trial  constraint  may  happily 
coincide  with  the  true  constraint  rQ  .  This  procedure  is  justified  by 
Theorem  1  if  the  hyperplane  given  by  zq'  and  that  given  by  kQdQ  are  in 
the  same  polyhedron. 


Remark  to  Stej 


A  more  detailed  procedure  to  identify  a  new  FOG  after 

C  X  <  (k  -  l)d  becomes  infeasible  is  as  follows: 

=  o  o 

If  d  >  2  ,  then  modify  the  added  constraint  C  X  <  k  d  -  1  for  trial 
o  ■  ■  o  o 

and  examine  whether  it  is  feasible  without  a  basis  exchange. 

If  it  is  feasible,  then  repeat  the  same  procedure  by  reducing  the  right- 
ham.  side  of  the  added  constraint  by  1  until  it  becomes  infeasible.  If  it 
is  infeasible,  then  a  new  G  H  F  can  be  identified  in  the  same  manner  as  in 
Remark  to  Step  2  and  Step  3  after  basis  exchanges.  The  trial  constraint  may 
happily  coincide  with  the  true  constraint  r^  . 

Thus  if  the  next  subpolyhedron  is  small,  we  may  "skip"  it  without 
Jumping  over  any  lattice  points  (Figure  4).  This  occurs  because  a  sub- 
polyhedron  cannot  contain  a  lattice  point  (from  Theorem  1)  if  it  does  not 


0 
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have  any  hyperplane  which  has  an  integer  value  of  the  objective.  In 

Figure  4,  a  subpolyhedron  *8  s^PPe<^  and  a  subpolyhedron 

V-VCV,V_  is  examined. 

J  J  D  / 
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4.  FINITENESS  PROOF 

Let  each  hyperplane  be  specified  by  Indices  K  ,  1  where  K  denotes 
each  subpolyhedron  and  1  denotes  each  hyperplane  in  the  same  subpolyhedron. 
(1  ■  0  for  K  such  that  the  K1"*1  polyhedron  does  not  contain  a  hyperplane 
given  by  a  multiple  of  dk  . ) 

Theorem  4.  (Theorem  of  Finiteness) 

The  algorithm  terminates  in  a  finite  number  of  steps. 

Proof ; 

The  number  of  i  (i.e.,  the  number  of  times  the  constraint  is  modified) 
for  each  K  is  clearly  finite  since  d  >  1  ,  and  K  is  also  finite 

l\  m 

because  the  number  of  the  extreme  points  is  finite  unless  the  linear  problem 
itself  is  unbounded.  For  each  i  ,  K  ,  the  number  of  basis  exchanges  is 
finite  because  a  hyperplane  is  n-1  dimensional.  After  identifying  the 
range  of  the  values  which  each  variable  may  assume,  the  enumeration  of  all 
permissible  combinations  of  the  integral  values  is  finite. 


b.  A  GRAPHIC  EXPLANATION  (FIGURE  5.) 


Step  1. 

We  obtain  Pq  ,  the  optimal  point  with  z  -  zq  .  Examine  whether 
Pq  is  a  lattice  point  or  not.  It  is  not  a  lattice  point  in  this 
case,  so  we  pass  to  Step  2. 

Step  2. 

Suppose  i  0  (mod.  d^)  for  all  i  . 

Step  3. 

Cut  the  convex  by  line  (1) 


C  X  <  k  d 
■  o  o 

where  dQ  is  the  greatest  common  divisor  of  all  the  variables 
belonging  to  FOG,  and  kQ  is  the  maximal  integer  such  that 
kQ  •  dQ  <  zq  .  Suppose  we  obtain  an  optimal  point  . 

Step  4. 

In  this  case,  it  is  feasible.  So  pass  to  Step  5. 


Step  5. 


By  exchanging  the  basis,  we  obtain  the  other  optimal  extreme 


point  Q 


(0) 

2 


Step  6. 

We  search  for  a  lattice  point  on  the  segment  with¬ 

out  success. 


Step  7. 


We  pass  to  Step  8  because  we  have  failed  in  obtaining  a  lattice 
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point. 

Step  8. 

Replace  line  (1)  by  line  (2) 

C  X  <  (k  -  l)d 
*o  o 

Solve  the  modified  problem  subject  to  (2)  without  the  basis 
exchange. 

Step  4. 

It  becomes  Infeasible  without  the  basis  exchange  because  we 
obtain  as  a  solution  the  intersecting  point  of  line  (2)  with 
line  (5). 

Step  9. 

We  identify  a  new  FOG  below  by  replacing  line  (2) 

by  line  (3):  C  X  <  k  d  -  3  to  obtain  a  new  G  after 

*=  o  o 

C  X  <  k  d  -  1  and  C  X  <  k  d  -  2  are  found  to  be  still 
"  o  o  *  o  o 

feasible. 

Step  2. 

Suppose  b^  =  0  (mod.  d^^)  for  all  i  again. 

Step  3. 

Modify  ro  into  r^  such  that 

C  X  <  k^ 

where  d1  is  the  greatest  common  divisor  of  all  the  variables 
elonging  to  the  current  set  of  FOG,  and  k^ 


is  the 
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maximal  integer  such  that  C  X  <  k^d^  is  in  the  subpolyhedron 

^0^1^2^3  ’  °^ta^n  as  a  s°luti°n  to  the  modified 

constraint. 


Step  4. 


It  remains  feasible  without  exchanging  the  basis. 


Step  5. 


We  obtain  the  other  intersecting  point  besides 

which  was  obtained  in  Step  3. 


Step  6. 


We  search  for  a  lattice  point  on  the  segment 


Step  7. 


A  lattice  point  is  not  obtained. 


Step  8. 


We  change  the  added  constraint  r^  to  r^ '  (line  (4)) 


C  X  <  (kx  -  l)dx  , 


and  we  obtain  as  a  solution  without  exchanging  the  basis. 


Step  4. 


The  solution  remains  feasible. 


Step  5. 


We  obtain  the  other  intersecting  point  7  by  introducing  the 
dual  degenerate  variable  into  the  basis 
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‘/.’VC 


Step  6. 

(2)  (2) 

We  search  for  a  lattice  point  on  the  segment  . 

Step  7. 

We  succeeded  in  obtaining  a  lattice  point  on  the  segment, 
successfully  terminating  the  algorithm. 


s*  ’V-'W  ' 


J 
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6.  NUMERICAL  EXAMPLE 


Maximize  z  *>  7x^ 

+ 

14x2 

+ 

10x3 

9x. 

+ 

5x~ 

+ 

3x. 

< 

34 

1 

2 

3 

4x, 

+ 

5x2 

+ 

6x^ 

< 

17 

1 

4x. 

+ 

10x„ 

+ 

7x_ 

< 

30 

1 

2 

3 

* 

Step  1. 

We  obtain  as  the  optimal  feasible  solution 

xi  ’  i*0  •  *2  "  2,6  *  x3  “  0* 


z  -  43.4 

The  solution  above  is  not  integral. 

Step  2. 

Add  a  constraint  for  trial 

7x.  +  14x„  +  10x_  <  43 
1  2  3  * 

partly  to  identify  the  variables  belonging  to  FOG. 

Step  4. 

The  solution  is  feasible. 

Step  9  and  Step  2. 

After  the  basis  exchange,  we  are  in  position  to  identify 
,  x2  and  all  three  slacks  as  the  variables  belonging 


to  G  D  F  . 
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d(1)  -  G.C.D.  (9,5,1)  -  1  , 
34  =  0  (mod.  1) 

,(3) 


d(2)  -  G.C.D.  (4,5,1) 
21  i  0  (mod.  1) 


1  , 


dw/  -  G.C.D.  (4,10,1)  -  1 
30  =  0  (mod.  1) 


Thus,  a  lattice  point  can  exist  because  the  slacks  In  G  make 
d^  *  1  for  all  1  in  this  case. 


Step  3. 

The  greatest  common  divisor  of  the  coefficients  of  and 

is  7  and  the  maximal  multiple  of  7  not  exceeding  43  is 
42  .  So  we  add  a  constraint 


7x^  +  14x2  +  10x^  <  42  , 


and  we  obtain  a  solution 


X1  “  3  * 


V  ■  — 

2  3 


and 


Step  4. 

The  solution  above  is  feasible. 

Step  5. 

By  introducing  one  of  the  dual  degenerate  variablep  into  the 


basis,  we  obtain 


3 


x 


1 


-  0 


and 


x 


3 


0 


We  successfully  obtained  an  optimal  Integer  solution. 


Identify  FOG  in 
a  subpolyhedron. 
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7.  A  METHOD  FOR  FINDING  LATTICE  POINTS 

A  problem  of  searching  for  lattice  points  in  Step  6  is  left  unsolved. 
One  method  is  to  repeat  the  same  procedure,  i.e.,  to  reduce  the  problem 
again  to  a  simpler  problem  to  find  a  feasible  solution  in  n-2  dimensions. 
At  this  time  d“l  is  uniformly  used  for  a  dummy  objective  function 
z  ■  £  x  in  all  subpolyhedrons;  we  do  not  have  :o  identify  G  H  F  anymore 

j  J 

because  the  value  of  the  objective  is  already  fixed. 

Another  method  is  to  use  a  branch  and  bound  algorithm  which  will  be 
presented  below.  Emphasis  is  placed  on  a  0-1  problem. 

As  stated  in  Theorem  2,  all  extreme  points  Q^,^,  ’  *  ‘  a  hyper¬ 

plane  H  can  be  obtained  by  repeating  the  introduction  of  variables  into 
the  basis  whose  shadow  prices  are  zero.  Let  ...,  x^\  •••* 

Xm+n)  denote  the  coordinate  of  (j  *  1 . N)  . 

If  there  is  any  lattice  point  on  H  whose  coordinate  is 

(x,,x_ . .  ...,  x  .  )  where  x,  should  be  an  integer  for  all  i  , 

1  l  1  nrt-n  1 

then 


*  min  (x|^)  <  xi  <  max  (x[^)  "  ui 


V  i 


If  [ ^ ”  [u  ]  <  ,  then  there  cannot  be  any  integer  x^  between 

and  . 

Rule  1.  (Criterion  of  Absence  of  a  Lattice  Point) 

If  [u  ]  <  hold  for  some  i  ,  then  H  In  disregarded. 

Rule  2.  (Criterion  of  Zero  Value  of  a  Variable) 

If  u ^  <  1  ,  then  x^^  ■  0  . 

Hereafter,  we  consider  only  0-1  case. 
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Rule  3.  (Criterion  of  Value  One  of  a  Variable) 

If  £  >  0  ,  then  x.  =  1  .  Rules  2  and  3  fix  some  variables 

X  j. 

on  H  . 


On  a  particular  hyperplane  H  ,  the  value  of  the  objective  function  is 

already  fixed.  Hence,  we  no  longer  care  about  optimizing  it.  The  order  in 

which  the  variables  are  assigned  values  is  unimportant,  contrary  to  the 

ordinary  branch  and  bound  algorithm  in  which  we  have  to  take  into  account 

the  variables  in  a  predetermined  order.  In  our  case,  the  only  criterion  that 

must  be  satisfied  is  \  c  x  =  k  d  .  Hence  we  will  take  the  assignment 

jeFDG  i  J 

of  the  variable's  value  in  the  simplest  order,  namely  in  the  ascending  order 
of  the  indices  of  the  variables. 

Hereafter,  we  assume  c^  >  0  for  all  j  .  (If  not,  we  can  obtain  it 
by  resetting  x^  =  1-x^  *  ^  F°r  t*ie  sa^e  simplicity  let  us  assume  the 

first  f  variables  are  already  fixed. 


An  Algorithm  of  Enumerations  by  Branch  Method 
n' 

r 

If 


(1) 


£  c  x.  *  k  d  for  n'  <  n  ,  then  we  obtain  an  integer 
j“l  j  J 

solution  by  assigning  value  zero  to  all  other  unfixed 
variables. 

(1-1)  If  the  solution  is  feasible,  we  have  obtained  an  integer 
feasible  solution. 

(1-2)  If  the  solution  iB  not  feasible,  then  disregard  this  branch 

and  search  another  branch, 
n' 

(2)  If  £  c.x.  <  k  d  for  n'  -  n  ,  then  go  back  to  the  branch 

j-1  i  J 

and  search  another  branch  for  a  larger  value  of  £  c  x  • 
n'  j  J  J 

(3)  If  £  c  x  >  k  d  for  n*  <  n  ,  then  go  back  to  the  branch 

j-1  J  3 

and  search  another  branch  for  a  smaller  value  of  )  C4X4  • 

J  J  i 
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A  Numerical  Example  of  Enumc  ..atirns 

Let  the  objective  function  be  9x^  +  I6X2  +  4x^  +  12x^  .  Let  us  assume 
x^  never  appears  as  a  basic  variable  ac  this  stage.  Hence  d  *  G.C.D. 
(16,4,12)  *  4  .  Suppose  k*5  is  appropriate  at  this  stage,  therefore 
z  -  5  x  4  -  20  .  Let  ^  -  (0,0.25,1,1)  ,  Q£  =  (0,0.875,0,0.5)  , 

Q-j  *  (0,1, 1,0)  ■  .  Thus,  0  <_  x^  <_  0  ,  0.25  £  X2  1  ,  0  ^  x^  ^  1  * 

0  _<  x^  _<  1  .  By  Rule  2,  x^  *  0  and  by  Rule  3  X2  =  1  .  The  only  unfixed 

variables  are  x^  and  x^  .  In  Figure  7,  a  figure  on  the  left  side  of  a 
slash  denotes  a  value  taken  on  by  a  variable  in  the  right  column  and  a 
figure  on  the  right  side  of  the  slash  denotes  a  value  of  the  objective 
function  taken  on  a  solution  associated  with  an  arc. 

We  start  with  Arc  1  to  which  a  value  16  is  associated  as  the  objective 
function.  16  is  smaller  than  20,  therefore  we  proceed  until  we  reach  or 
exceed  20.  We  consider  Arc  2  where  we  reach  an  end  of  a  branch,  but  its 
value  16  is  too  small,  so  we  search  another  branch  for  a  bigger  value. 
(Procedure  2.)  We  consider  Arc  3  but  its  value  28  is  too  large,  therefore 
we  go  back.  (Procedure  3.)  We  consider  Arc  4  where  its  value  is  exactly 
20.  Thus,  we  assign  x^  a  zero  value.  Hence,  automatically  we  move  to 
Arc  5  by  disregarding  Arc  6.  (Proceoure  1.)  We  examine  whether  Arc  5  is 
feasible  or  not.  Arc  5  happens  to  coincide  with  ,  therefore  it  is 
feasible.  (Procedure  1-2.) 

Another  Numerical  Example 

Let  us  consider  the  same  objective  function  and  a  different  k=6  . 
Hence,  the  criterion  is  9x^  +  I6X2  +  4x^  +  12x^  =  24  .  Let 
Qj_  -  (0,1, 0,2/3)  ,  Q2  =  (0,1, 1,1/3)  and  Q3  =  (0,1,1/2,5/12)  =  Q4  . 
u^  =  2/3  [u4]  =  0  <  =  1/3  .  Hence,  there  cannot  be  any  integer  feasible 


solution.  Thus,  we  must  disregard  this  hyperplane.  (Rule  1.) 
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FIGURE  7:  A  TREE 


fixed 


fixed 


(0,1) 


{0,1} 
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The  algorithm  presented  above  is  a  branching  algorithm;  but  it  is  a 
bounding  algorithm  in  a  somewhat  different  sense  than  usually  conveyed  by  the 
term  "branch  and  bound."  For  each  H  a  value  of  the  objective  function  is 
already  fixed,  so  we  do  not  need  to  evaluate  it  anymore.  Instead,  we  give 
an  upper  and  a  lower  bound  to  each  variable. 

A  criterion  t  nction  given  here  is  the  criterion  to  see  whether  an  arc 
is  feasible  or  not,  as  opposed  to  the  ordinary  branch  and  bound  algorithm 
where  the  criterion  is:  which  variable  should  be  considered  next. 

No  evaluation  after  setting  a  hyperplane  H  is  made  in  the  algorithm 
presented  above. 
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8.  DISCUSSIONS 

Another  Starting  Procedure 

We  have  started  with  an  optimal  continuous  solution.  But  we  can  also 
start  with  a  feasible  integer  solution  if  it  is  already  given.  In  the  latter 
case,  we  always  have  a  feasible  integer  solution  and  we  proceed  to  a  better 
integer  solution  until  an  optimal  continuous  solution  is  reached  or  exceeded. 
The  termination  criterion  is  the  same  as  in  the  former  case. 


The  Possible  Application  of  the  Algorithm 

Let  us  assume  that  the  matrix  A  (in  A  X  “  B)  is  m  x  (n-Hn)  (with 

m  slack  variables).  Without  degeneracy,  m  basic  variables  are  associated 

at  each  0^  ,  and  one  extreme  point  is  already  identified  on  H  .  Hence, 

(n-1)  nonbasic  variables  assume  zero  shadow  prices  at  .  Thus, 

(m  +  (n-1))  variables  belong  to  at  .  This  leaves  out  one  variable. 

In  particular  if  a  slack  variable  is  this  exception,  then  F  0  •  F  .  But 

consider  the  problem  of  m  constraints  D, (x, . x  )  >  0  (i  ■  1,  ...»  m)  , 

i  1  n  ■ 

of  which  at  least  m'  constraints  are  required  to  be  satisfied  with  Integral 


Xj  for  all  j  .  The  problem  is  represented  as  follows: 


Di  "  Vi  *  0  (i  -  1 . m) 

A^  -  m  -  ra* 

A^  “  (0,  1} 


where  for  each  i  .  In  this  problem  in  addition  to  (m+n) 

variables,  there  are  ra  A-variables  whose  coefficients  are  all  zero  in  the 

objective  function.  The  cardinal  number  of  G.  ■  m+n  -  1  >  the  cardinal 

i 

number  of  F  D  G^ 


>  max  (0,  n-m  -  1) 


0 


while  the  total  number  of  the  variables  equal  2  m  +  n  .  Thus  we  can  see  if 
n  i ;  not  too  large  for  m  ,  the  set  F  H  is  a  relatively  small  subset 
of  the  set  of  all  the  variables. 

When  G  consists  mostly  of  slack  variables  and  A^'s  •  d  may  be  greater 
than  one.  Thus,  a  cutting  hyperplane  method  may  work  efficientl>.  On  the 
other  hand,  many  d^'s  are  equal  to  one  because  coefficients  of  slack 
variables  are  one  in  constraints. 

When  G  consists  mainly  of  nonslack  variables,  then  the  d^'s  may 
be  greater  than  one  and  S^ep  2  may  work  efficiently  in  excluding  from 
consideration  subpolyhedrons  which  cannot  contain  a  lattice  point.  On  the 
other  hand,  the  cutting  hyperpL.ane  method  may  not  work  so  well  in  this  case 
because  d  may  be  close  to  cne. 
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